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Abstract: A modification to the contransverse mass (Mct) technique for measuring the 
masses of pair-produced semi-invisibly decaying heavy particles is proposed in which Mct 
is corrected for non-zero boosts of the centre-of-momentum (CoM) frame of the heavy 
states in the laboratory transverse plane. Lack of knowledge of the mass of the CoM frame 
prevents exact correction for this boost, however it is shown that a conservative correction 
can nevertheless be derived which always generates an Mct value which is less than or 
equal to the true value of Mct i n the CoM frame. The new technique is demonstrated 
with case studies of mass measurement with fully leptonic tt events and with SUSY events 
possessing a similar final state. 
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1. Introduction 

Techniques for measuring the masses of pair-produced particles decaying semi-invisibly 
through short decay chains at hadron colliders have attracted considerable interest. The 
principle motivation for the development of such techniques is the measurement of the 
masses of supersymmetric particles ('sparticles') at the Large Hadron Collider |l], ^, |3|, f|, |5|, 
|6|, [^, |8|, |9|, [l(], 12 , 12 , 14 , 15 , [l6| , however they may be applied more widely to measure 



the mass of the top quark at the Tevatron |17| or LHC Pq |, or to identify fully leptonic 
WW events 0. 

Recently || a straightforward new variable, the 'contrans verse mass' (Mct), was pro- 
posed which enables the measurement of a simple analytical combination of the masses of 
the pair-produced heavy states 5i (i = 1,2) and their invisible decay products Oj. The 
contransverse mass is defined by 

M CT (v 1 ,v 2 ) = [E T (vi) + E T {v 2 )] 2 - [pT-(fi) -Pt(^2)] 2 
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= m 2 ( Vl ) + m 2 (v 2 ) + 2[E T {v 1 )E T {v 2 ) + p T («i) • Pt(^)], (1.1) 

where Vi are the visible products of each decay chain, p-rC^i) is the tranverse momentum 
vector of Vi and 

Er(vi) = yjrtivj+mtivi). (1.2) 

It can be shown |J that Mct is in general bounded from above by a quantity dependent 
upon the masses m(5) and m(a). If m(vi) = m{v2) = m(u) then the distribution of event 
Mct values possesses an end-point at: 

M^\ m \»)\ = ^ + m2{S) -f {a) . (1.3) 
m(oj m(o) 

Consequently a measurement of the gradient and intercept of the linear function describing 
the dependence of M^ x on m 2 (v) allows both m(5) and m(a) to be measured indepen- 
dently. 

Despite the simplicity and ease-of-use of the contransverse mass technique, it suffers 
from two principle draw-backs p|. The first is that Mct is not invariant under Lorentz 
boosts of the 6162 centre-of-momentum (CoM) frame in the laboratory transverse plane. 
Consequently if the <5ic>2 system recoils in the transverse plane against upstream object(s) 
such as ISR jets then the value of Mct calculated in the laboratory frame is not in general 
equal to that calculated in the 6162 CoM frame. Mct values can be generated which are 
greater than M™|; x and as a result the Mct end-point can be smeared (see e.g. Figure 2 in 
Ref. 0). The second draw-back is apparent when attempting to measure m{5) and m(a) 
independently in events with non-zero visible masses m(vj) using Eqn. ( |1.3|) . The require- 
ment m{v\) = m{v2) can significantly reduce the available event statistics and require the 
accumulation of very large integrated luminosity, even for channels with relatively large 
a.BR. This problem is illustrated clearly in Figure 3 in Ref. Q. This paper will seek to 
address these two problems and demonstrate the utility of the Mct technique through two 
case-studies. In the process we shall identify a further problem with using Eqn. (1.3) to 
measure masses independently, but develop an alternative strategy for two-step sequential 
two-body decay chains combining Mct end-point measurements with conventional invari- 
ant mass end-point constraints. We shall also investigate the use of the transverse boost 
dependence of M™|; x to measure masses independently, but find that this technique suffers 
from similar problems. 

The structure of the paper is as follows. Section ^ will study the transformation 
properties of Mct under contra- linear and co-linear Lorentz boosts of Si, leading to the 
development of a procedure for correcting Mct for co-linear boosts. Section || will discuss 
the shape of the resulting Mct distributions. Section || will propose a new method for 
maximising the available event statistics when measuring m{5) and m(a) independently 
with Eqn. (1.3) by removing the m{v\) = m{v2) requirement. This section will also develop 
a technique through which m{5) and m(a) can in principle be measured independently by 
using the transverse boost dependence of M™|; x . Section [| will investigate these techniques 
with LHC case studies measuring the masses of the top quark, W and neutrino with fully- 
leptonic tt events, and the masses of SUSY particles decaying to a similar final state. 
Section || will conclude. 
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2. Transformation properties of M C t 



2.1 Equal magnitude contra-linear boosts 

It is instructive to consider first the transformation properties of Mct under contra-linear 
equal magnitude boosts, in which b\ and b 2 move in opposite directions with equal mo- 
mentum. Mqt is derived from the quantity Mc given by 

M ( Vl ,v 2 ) = [E( Vl ) + E(v 2 )} 2 - [p( Vl ) - P (v 2 )} 2 

= m 2 ( Vl ) + m 2 (v 2 ) + 2[E{ Vl )E{v 2 ) + p(v a ) • p(v 2 )}, (2.1) 

and this was shown in Ref. || to be invariant under such boosts. By contrast Mqt is n °t 
in general invariant under such boosts, however the position of the Mct end-point, Mq^ x , 
is. The reason for this can be understood by observing that one can rewrite Eqn. fl2.1|) in 
the following form in the <5i<52 CoM frame: 

M c {v x ,v 2 ) = m 2 {v l ) + m 2 {v 2 ) + 2[E T {v l )E T {v 2 )coshT 1 r,{v i ) + p T {v l )-p T {v 2 )}. (2.2) 

Comparing with Eqn. Ql . 1| ) and noting that cosh £77(1^) > 1 this shows that Mct < Mc, 
with equality when £77(1^) = 0. Now Mc, like Mct, is bounded from above by M™|; x and 
so one finds finally that Mct < Mc < M™^ x . It is interesting to note additionally that 
Mc equals M^ x when v\ and v 2 are co-linear in the Si rest frames and hence the necessary 
and sufficient criteria for Mct = Mqj? are that At](vi,v 2 ) = and A(f>{v-\_,v 2 ) = in the 
5i rest frames while £r/(-Uj) = in the <5i<!>2 CoM frame. 

A similar argument applies when the transverse mass Mt undergoes co-linear equal 
magnitude boosts, and the result is similar, namely that Mt is not invariant under arbitrary 
transverse boosts but nevertheless possesses a boost-invariant end-point. It is interesting 
to note that in this case the necessary and sufficient conditions for Mt = M™ ax are in 
some sense the complement of those in the Mct case: here £77(1^) = and Acj)(vi,v 2 ) = it 
in the rest frame(s) of the parent particle(s) while Ar](vi,v 2 ) = in the event CoM frame. 
Of course if v\ and v 2 are the sole products of the decay of the same parent then the first 
and second criteria are generally satisfied through conservation of momentum. 

2.2 Equal magnitude co-linear boosts 

The contransverse mass is invariant by construction under co-linear equal magnitude boosts 
of 5\ and 8 2 in the beam (z) direction, by virtue of its dependence purely on transverse 
quantities. Consequently Mqj? is similarly invariant. In the presence of co-linear equal 
magnitude boosts in the CoM transverse plane however, equivalent to a single global trans- 
verse boost, the values of both Mct and M™|; x can depend on the magnitude and direction 
of the boost. For instance, when m{v\) = m(v 2 ) = m(v) Eqn. (|1.3[) becomes: 



M^[m 2 (v), Pb } = 2(rp + EoVT+^), (2.3) 

where r = pb/2m(5), pb is the net transverse momentum of upstream objects (ISR jets 
etc.) generating the boost, 

7-1 _ m 2 (5) -m 2 (a) + m 2 (v) 

Eo = 7, — 77\ 1 \ ZA ) 

2m(d) 
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and po = Eq — m 2 (v). 

With sufficient statistics one might hope to use Eqn. (2.3) to measure m{5) and m(a) 
separately by measuring Mqj? as a function of pt,. This possibility is considered further in 
Section With limited statistics however it would be useful to be able to transform Mct 
such that its value always lies below the Mqt end-point given by Eqn. ([D^). In this case 
one sacrifices the ability to measure m(5) and m(a) independently using the p^ dependence 
of Mqj? in order to maximise statistics near the pb = Mct end-point, while limiting 
smearing beyond the end-point due to integration over p^. 

The approach we shall take will involve boosting the four-momenta of the visible decay 
products Vi back into the 5\5 2 CoM frame with boost factor /3, prior to calculating Mct- 
If we know neither the sign nor the magnitude of f3 then the minimum value of Mqt that 
we can obtain by varying the assumed value of f3 is the one-dimensional analogue of Mct 
given by Mc y defined by 

M% y { Vl ,v 2 ) = [E y ( Vl ) + E y (v 2 )] 2 - [ Py ( Vl ) - p y (v 2 )} 2 

= m 2 ( Vl ) + m 2 (v 2 ) + 2[E y ( y v 1 )E y (v 2 ) + Pyiv^py^)}, (2.5) 

where E 2 (vi) = p 2 {vi)+m 2 (vi) and we have assumed that the boost lies in the ±x direction. 
Mc y is invariant under arbitrary boosts in the x — z plane for the same reason that Mct 
is invariant under boosts in the z direction. It represents a conservative lower bound on 
the value Mct(CoM) of Mct measured in the &\b~ 2 CoM frame. The criterion for Mct to 
equal Mc y in any given frame, and hence for Mct to be minimised in that frame, is that 
A x = 0, where A x is defined by 

A x = Px{vi)E y {v 2 ) + p x (v 2 )E y (vi). (2.6) 

As an aside, it is interesting to consider at this point the possibility of using Mq defined 



by Eqn. (2.1), rather than Mct, and attempting to perform a correction for longitudinal 
boosts along the beam direction. In this case we know neither the sign nor the magnitude 
of the z-boost and so by the above argument the appropriate quantity to use is the two- 
dimensional analogue of Mc, which is just Mct as we have been using already. The 
criterion for Mc to equal Mct is by analogy with Eqn. (|2.6[) 



P z { Vi )Et{v2) + Pz(v 2 )E T ( Vl ) = 0, (2.7) 



which is equivalent to setting Y,r](vi) = as required by Eqn. ( |2.2| ). 

Now in fact we do know the sign of the required boost, because we know the direction 
in the transverse plane of the upstream momentum. Defining this direction to be the +x 
direction we need to boost v\ and v 2 in this same direction (i.e. use /3 > 0) to correct 
for the original boost of the CoM frame, which must have been in the — x direction. 
Such a +x boost monotonically increases p x {vi) and hence monotonically increases the 
transformed value of A x towards +oo as f3 — > +1. Let us define now A x n a ^ and M C t riab) 
to be the values of A x and Mct measured in the lab frame and A' x and M C t(cott) to be 
the equivalent values obtained after boosting v\ and v 2 . If A x (\ &h ^ > then M C T(corr) 
increases monotonically from M C T(\ab) towards +oo as f3 — > +1 (see Figures [j] and |2|). 
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Figure 1: Schematic diagram showing the dependence of Mc T / COIT \ on the boost factor /3 used 
in the boost correction. Cases with A x n a ^ > and A x ^ a ^ < are shown. When /3 > and 
Ar(iab) > the minimum value of McT(corr) occurs when (3 = and hence Mct{coh) = -^CT(iab)- 
If p > and ^4 x (i a b) < the minimum value of AfcT(corr) is Mc y - 

Consequently in this case the least conservative lower bound on AfcT(CoM) we can obtain 
is M(jx{cow) = ^CT(iab) • If on the other hand ^4 x (i a b) is negative then as /3 and A' x increase 
-^CT(corr) nrs t decreases from McT(iab) towards its minimum value of Mc y (at A' x = 0) 
before increasing again towards +oo (see Figures | and |). In this case, without further 
information, the best we can do is set -McT(corr) = ^Cy 

Fortunately however we have not yet exhausted the possibilities for boost correction. 
Observe first that when boosting v\ and V2 the boost factor is given by f3 = pb/Egf where 
Egf 1 is the assumed value of the energy Ess of the 5\52 CoM frame in the lab frame. 
Consequently increasing the value of (3 is equivalent to decreasing Egf, and vice versa. 
Hence if Mct(cott) increases (decreases) monotonically with increasing Egg and we set 
Egf to a value less than (greater than) Ess, the value of M(j T r COII \ we obtain provides a 
conservative lower bound on M CT ( CoM ^ . Now if A x /i ah j > then M CT ( C orr) always increases 
with increasing (3 (see Figure H) and hence it decreases with increasing Egf (see Figure || 
- upper curve). In this case we should set Egg to the upper bound on Ess, boost v\ and 
V2, and obtain a conservative lower bound on Mct^oM) from the value of M^^( corr ) in this 
frame. If A x ^ ah ^ < the situation is more complicated (see Figure || - lower curve). In this 
case, if A' x < after boosting with Egf set to both the upper and lower bounds on Ess 
then the least conservative lower bound on M CT ^ CoM ^ is given by the value of M CT ( CO rr) 
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Figure 2: Schematic diagram showing the dependence of Mct(cott) on A' x . If Ar(iat>) < (left- 
hand figure) then as A' x increases from A x ^ a ^, AfcT(corr) passes through its minimum at A' x = 0, 
while if ^(lab) > (right-hand figure) it does not. 

with Egg set to the lower bound on Egg. Conversely if A' x > in both these cases then the 
least conservative lower bound on Mct(c oM ) is given by the value of M CT ^ COIV ^ with Egg 
set to the upper bound on Ess- If A' x > with Egg set to the lower bound on Egg but 
A' x < with Egf set to the upper bound on Ess then Mct(CoM) could be as low as Mc y 
and so this should be used as the least conservative lower bound on Mct(CoM) ■ 

In fact we can indeed obtain both upper and lower bounds on £"5,5. An upper bound 
is provided by the proton-proton centre of mass energy E cm while the total visible energy 
E of the decay products provides a lower bound. This latter quantity is calculated by 
summing the energies of the visible decay products with the net transverse momentum of 
the invisible decay products 1 given by EJp lss . EJp lss equals the total energy of the invisible 
decay products only when these are massless, co-linear, and moving in the transverse plane, 
and so in general E < Ess- Below we shall denote values of A' x obtained with Egg set to 
E cm or E as respectively ^ (lo) and A' x( ^ . 

Let us now summarise the procedure we have developed for correcting Mqt for the 
effects of co- linear equal magnitude boosts of S± and 82 in the transverse plane 2 . First 
calculate A x ^ ah ^ and A',^ using Eqn. (p.6| ), the latter by boosting v\ and V2 with /3 = 
Pb/E cm . If A^iab) > or A' x Q o -. > then one should set M CT ^ corr ^ to the boosted value of 
Mqt obtained with /3 = pb/E cm . If neither of these criteria are satisfied then one should 
next evaluate E and boost v\ and V2 with f5 = pb/E. Now evaluate ^,(hi) usin g Eqn. ([O]). 

If a lower bound m\ (a) on the masses of the individual invisible decay products can be assumed then 
conservatively one can use y 7 (_E™ 1SS ) 2 + 4mj^(a) in E instead of _E™ 1SS to obtain an improved bound on Ess- 
2 f77, CH — h and ROOT code implementing this boost-correction procedure can be downloaded from 
http : //projects . hepf orge . org/mctlib. 
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E ss <*t (Arb. Unils) 



Figure 3: Schematic diagram showing the dependence of M CT ( COIT ) on Egf, the estimated value 
of Ess used in the boost. Cases with Ar(iab) > and A x n & \,-\ < are shown. When Awiat,) > or 
A' x > M C T(\ab) decreases with increasing Egf. If however A x n a ^ < and A' x < then McT{\ab) 
increases with increasing Egf. 



If A',^. < then one should set Mct(coti) to the value of Mqt in this boosted frame. If 
however A',^ > then one should set Mqt(coh) = Mcy An example of the effect of this 
boost correction procedure is shown in Fig. |] for the SUSY events considered in Ref. || . 



3. The shape of the Mqt distribution 



The differences between My and Mqt identified in Section 2.1 affect the shapes of the 
distributions of these quantities. As is well-known, Mt possesses a Jacobian peak at 
Mt = M™ ax when v% and V2 are the sole products of the decay of the same parent. 
Physically this peak arises because near the end-point all kinematic configurations with 
different rj(vi) generate the same value of Mt, in other words Mt becomes independent of 
the kinematics of v\ and V2- 

Turning now to Mqt, let us consider first the special case where the 6182 system is not 
boosted in the laboratory transverse plane and no boost correction is applied. When Mqt 
is calculated for the visible decay products of the 8% 82 system the extra degrees of freedom 
resulting from the independent motion of v\ and V2 generate a significantly different shape 
of distribution. Near the end-point at M^ x only events in which both v\ and V2 move 
in the transverse plane can contribute to the distribution. The small probability of this 
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Figure 4: Mct distributions of SUSY events containing at least two jets for the SPSla benchmark 
SUSY model considered in Ref. ||. In the left-hand figure the open histogram shows the Mct 
distribution with no pb cut and no correction applied. This diagram should be compared with 
Figure 2 of Ref. [p| showing the evolution of the uncorrected Mct end-point as a function of the 
cut onpb. The light (yellow) histogram shows the Mct distribution of the same events after the co- 
linear boost correction described in the text has been applied to the Mct values, assuming that all 
additional jets contribute to the upstream momentum pb- The medium (cyan) histogram shows the 
same distribution for quqii pair-production events. In the right-hand figure the latter distribution 
is plotted on an expanded scale. The open histogram shows the parton-level Mct distribution for 
the same events. The end-point from q^qu pair-production is expected at 531 GeV (denoted by a 
vertical dashed line in both figures). 



configuration (because v\ and V2 are uncorrelated) cancels the large probability generated 
by the Jacobian transformation, resulting in an end-point which tends asymptotically in 
the absence of boosts to 

P(M CT ) dM CT = A^J(M^ X ) 2 - M CT dM CT , (3.1) 

where A is a constant. Typical Mqt distributions in the absence of boosts, displaying this 
end-point, are shown in Figure [|. 

Despite this cancellation of the peak at the Mqt distribution can still possess 

a Jacobian peak. The peak occurs however at the lower limit of the distribution, where 
Mct = m(v\) + m(v2) = M™™ (see Eqn. ( |1.1[ )). The distribution tends asymptotically in 
the absence of boosts to 

Mct 

P(Mct) dM CT = B = dM CT , (3.2) 

V m ct - mpy 

where B is a constant. This is very similar to the functional form of the My Jacobian 
peak, although in this case it is reversed such that the distribution is real above the peak 
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Figure 5: Typical Mct distributions in the absence of boosts for massless visible particles (full 
histogram) and massive visible particles of total mass 0.13M™|; X (dashed histogram). The x-axis 
has been scaled such that the end-point lies at Mct/Mqj? = 1. 



rather than below it. Physically the peak occurs because when Mqt ~ M™™ the value 
of Mct becomes independent of the kinematics of v\ and V2, because the m(vi) terms in 
Eqn. (IT) dominate. Hence all kinematic configurations generate similar values of Mct- 
An example of such a peak can be seen in the dashed histogram in Figure ||. Note that 
when Mc?p = the numerator and denominator in Eqn. ([^) cancel leaving a uniform 
distribution (see e.g. full histogram in Figure |5|). 

Consider now the general case where the 8182 system has been boosted in the laboratory 
transverse plane, and the boost correction procedure discussed in Section |2.2| has been 
applied. In this case there can be a further enhancement of the population of events at 
M C t = Mgj?. If A,.( lab ) < and A',^ > then the boost-corrected value of M C t is 
given by Mc y from Eqn. (|2.5| ). If the tranverse momenta of the two visible particles under 
consideration are bisected by the boost direction however, then Mc y can take very small 
values, even if the transverse momenta of the V( particles are relatively large. This effect 
is particularly striking when m(vi) = m(v%) = M™™ = 0, in which case it is straight- 
forward to see from Eqn. ( |2.5[ ) that Mc y = 0. We shall refer below to the resulting peak 
at M C t = M^ as the 'M C t = M Cy peak. 

Specific examples of the peaks and end-points discussed in this section can be seen in 
Section [| below. 
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4. Measuring m(5) and m(a) independently 



The dependence of M^f? on m ( v i)i noted in Ref. M, and on pt>, noted in Section 2.2 
provides potential techniques for measuring m{5) and m(a) independently. In this section 
we develop these techniques in more detail. 

4.1 Using the m(vi) dependence of M™|; x 

As remarked in Section |], when attempting to measure m{5) and m(a) independently 
using Eqn. ( |1.3| ) the requirement m{v{) = m(v 2 ) = m(v) reduces significantly the event 
selection efficiency. One can consider ameliorating this problem by removing this mass 
equality requirement and considering the dependence of the resulting Mqt end-point on 
both m 2 {v\) and m 2 (v 2 ). This is given by 

(M^# x [m 2 ( Vl ) , m 2 (v 2 )] ) 2 = m 2 ( Vl ) + m 2 (v 2 ) + 2 ( E ( Vl ) E (v 2 ) 

+ ^[£ 2 M - m 2 { Vl )] [E 2 {v 2 ) - m 2 (v 2 )]), (4.1) 

where 

m 2 {5)-m 2 (a) + m 2 (v i ) 

E (vi) = ^ . (4.2) 

In this case all events passing background rejection cuts are used, however the implicit 
requirement of binning in both m 2 {v\) and m 2 (v 2 ) to measure M™^[m 2 {v\) : m 2 (v 2 )] still 
limits the available event statistics in each end-point measurement (modulo the symmetry 
under interchange of v\ and v 2 of M^^[m 2 (v i), m 2 (v 2 )]). As an aside, Eqn. ( f4.ip provides 
a link between the 'stransverse mass' Mt 2 (x) 0] an d Mct- This is discussed in more detail 
in Appendix A. 

An alternative to using Eqn. ( |4.1| ) for measuring m{5) and m(a) independently involves 
observing that M^^\m 2 {v)\ is linearly dependent on m 2 (v) in Eqn. ( |l.3|) and hence that 

M c T[<axl = max (M^[m 2 { Vl )\, M^[m 2 (v 2 )]) , (4.3) 

where, 

m max = max(m(ui),m(u 2 )) • (4.4) 
Then we can make use of the following inequality: 

M^m 2 ^),™ 2 ^)] < M^fVLj, (4-5) 

to find that 

M CT (vi,v 2 ) < M^m 2 ^),™ 2 ^)] < MgTKj. (4.6) 

Consequently if the two-dimensional distribution of event Mct(vi, v 2 ) values versus event 
m max values is plotted, for all events passing background rejection cuts, the distribution 
will display an Mqt end-point dependence on given by: 



M maxr 2 i _ <ax . m 2 (5)-m 2 (a) 

M C T [m max \ - ^ + m{§) . (4.7) 
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Hence m{5) and m{a) may be obtained by measuring the gradient and intercept of the end- 
point dependence on m max in a similar manner to the existing technique using Eqn. (|1.3|). 
Now however all events passing the background rejection cuts can be used rather than just 
a small subset. 

It should be noted here that although this technique is sound from a theoretical point- 
of-view, the uneven distribution of events in the Mct(v\,V2) versus m max plane can cause 
difficulty when attempting to use it in practice. This is discussed further in Section [|. 

4.2 Using the pb dependence of Mqj? 



In Section ^ we developed a procedure for correcting Mct{v\,V2) such that it is always 
bounded from above by the expression for M™^ x obtained when the upstream momentum 
Pb = 0. Given sufficient statistics however an alternative procedure would involve binning 
the non-boost-corrected value of Mct(vi, v 2 ) in Pb and measuring m(S) and m(a) indepen- 
dently from the dependence of Mqj? on p&. When m(vi) = m(v2) = m(v) this dependence 
is given by Eqn. (^j^), however we can also obtain a general expression valid even when 
Tn(v{) ^ m(v2), which is 

(MSr[m 2 (^i),m 2 (^),p 6 ]) 2 = (M^[m\v 1 ),m 2 (v 2 )]y+ 

4r 2 (EolyJEofa) + Po(vi)p (v 2 ) + ^1 + r 2 [p (vi)E (v 2 ) + PoM^oM] ) , (4.8) 
where M^j^[m 2 (vi) , m 2 (v 2 )] is obtained from Eqn. (|4.1| ), r = pb/2m(5), Eo(vi) is given by 



Eqn. ([h^) andpo(vi) = y/ E'^vi) — m 2 {vi). In principle the additional dependence of M™^ x 
on the angle between the pb vector and the net momentum of v\ and Vi in the transverse 
plane could also be exploited, although this is not considered further here. 

The advantage of using the pb dependence of M™^ x rather than the m(vi) dependence 
discussed in Section |4.1| is that m(5) and m{a) can in principle be measured independently 
even when m{y\) = m(v2) = m(v) ~ 0, for instance when the Vi particles are jets or 
leptons. This avoids potential combinatorial problems inherent in the latter technique. In 
this special case the expression for Mqj? becomes 

M^[0, 0,p b ] = M^[0, 0, 0] (r + s/l + r 2 ) . (4.9) 

A toy Monte Carlo example Mct{vi,V2) versus pb distribution for massless Vi particles is 



shown in Figure o together with the theoretical bound from Eqn. 4.£. 



As with the use of the m(vi) dependence of M™|; x in Section 44 the above technique is 
sound from a theoretical perspective, however as we shall see in Section [| the uneven dis- 
tribution of events in the Mct( v i, ^2) versus pb plane can cause difficulty when attempting 
to use it in practice. 

5. Mass measurement for two-step decay chains 
5.1 Generic strategy 

We shall now investigate the use of the boost-corrected contransverse mass discussed in 
Section |2.2| (hereafter referred to simply as Mct) to measure the masses of pair-produced 
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Figure 6: Example Mct(vi, V2) versus pb distribution generated with a toy Monte Carlo showing 
the dependence of M™|; x on pb given by Eqn. (fb|). The simulated events consist of pair-produced 
5 particles of mass 200 GeV, each decaying into a massless visible particle and an invisible particle 
a of mass 100 GeV. pb is evenly distributed in the range 0-300 GeV. The Mct(w 1 ,w 2 ) axis has 
been normalised to M™|; x obtained when pb = 0. 



heavy particles decaying via symmetric two-step sequential two-body decay chains. As 
discussed in Ref. H] mass measurement with such chains is non-trivial because they are 



too short to solve fully for the masses using invariant mass end-point techniques [20, 21 



22, [23|, 24] or the 'mass relation' method ]25(| . The decay chains considered can be written 
in the form: 

5 -»• Pfi -»• PQa, (5.1) 

where 6, (3 and a are generic massive particles, P and Q are generic visible particles (here 
assumed massless) and a is invisible. The two chains present in each event can be seen in 
diagrammatic form in Figure [?], where particles appearing in the second decay chain are 
denoted with primed labels. We assume in the following discussion that particles labeled 
with the same letter possess the same mass. 

With each event we can construct one pair of invariant mass observables and three 
contransverse mass observables from the momenta of the four observed particles P, Q, P' 
and Q' . These observables are: 

• m(P", Q^'')' the invariant masses of the visible products of the two decay chains 

• Mqt(P, P')- Mqt constructed from the momenta of P and P' 
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• Mct{QiQ')'- Mqt constructed from the momenta of Q and Q' 

• Mct{[PQ], [P'Q']) : Mct constructed from the momenta of the aggregate products 
of each chain [PQ] and [P'Q'\. 

These observables possess kinematic end-points whose positions are functions of the masses 
m(5), m(/3) and m(a). The end-point positions are respectively 3 : 



MW/) = ^,i,, (5.3) 
A^(Q,9-) = m2( %f (a) ^fe3, (5.4) 
M^([P«], [P'QD = m '% m ' W + m(S> ( m2 %™ 2(a) ) . (5.5) 



where the final relationship is obtained from Eqn. ( |1.3| ) with m(v) = m ma,x (P,Q). In 
addition the two-dimensional distribution of events in the Mct([PQ], [P'Q']) versus 



max 



plane discussed in Section 4.1 can be constructed, providing additional mass constraints 
via Eqn. (fl.Tj). The two-dimensional distributions of events in the (non-boost-corrected) 
Mct(P, P') versus pb and Mct(Q, Q') versus pb planes discussed in Section can also in 
principle be used. 

Using Eqns. ( |5.3| ) and ( |5.4| ) the mass of the parent particle 5 can be calculated 

from: 



kfk 



m{6) = k^we- (5 - 6) 

The masses of j3 and a can then be obtained by simple substitution into Eqns. 
and (|5.4j) respectively. The constraints on the masses provided by Eqn. (p~7|), Eqn. (I 



and/or Eqn. ( |4.8D may be more difficult to exploit, as shall be discussed in Section but 
they can be used as a closure test for the masses measured using the other constraints. 

In typical SUSY decay chains we consider cases where the particles P and Q can be 
either quarks/jets (q) or leptons (£), leading to the following possible final-state config- 
urations: {P = q,Q = q}, {P = q,Q = £} and {P = £,Q = £}. The second of these 
configurations is particularly favourable from an experimental point-of-view because in 
this case there is no ambiguity in assigning particles to steps in the decay chains when 
constructing Mqt(P, P') and Mct(Q,Q')- In the case-studies presented below we shall 
therefore focus on events with the final state P = q and Q = £. 



When dealing with single-step three-body decay chains in which two visible particles and one invis- 



ible particle are produced in each decay Eqns. ( 5.2 )— ( 5.S ) are replaced by m max (P, Q) = m{8) ~ m(ct) 



M^([PQ],[P'Q']) = 2[m(6)-m(a)] and M^(P,P') = M$^(Q,Q') = [m 2 (S) -m 2 (a)]/m(5), although 
in the last two cases the distributions are strongly phase-space suppressed near the endpoints. 
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Figure 7: Diagrammatic view of the decay chain described in the text. The mass measurement 
technique described in the text is independent of the boost given to the system of interest by 
upstream decays or ISR (denoted by circles). 

5.2 Benchmarking on top events 

The mass measurement technique proposed above can be tested with ti events in which 
both top quarks decay to leptons via leptonically decaying W's through the chain: 

t -> bW -> blv. (5.7) 

Such events contain two symmetric two-step sequential two-body decay chains, with invis- 
ible particles being produced at the end of each chain. They therefore provide a suitable 
testbed for our mass measurement technique, with 5 = t, (3 = W, a = v, P = b and 
Q = I. The main notable difference between these events and SUSY events is that the 
invisible particles are in reality approximately massless, however in our analysis we shall 
not make this assumption. This approach has been used previously to study alternative 
SUSY particle mass measurement techniques 14]. 



In order to evaluate the observables discussed in Section 5.1, we generated with MC@NL0 



3.3 [26, G% an inclusive sample of ^/s = 14 TeV LHC tt events with an input top mass of 



172.5 GeV. The events were passed through the parameterised detector simulation ACERDET 



1 28] which was modified to reproduce the resolutions for leptons and jets given in [29]. For 
the tagging of 6-jets, an efficiency of 60% was assumed, for a light jet rejection of 100. 
The total generated sample was 2.2 M events, corresponding to an integrated luminosity 
of approximately 3 fb . 
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Events were selected with the following requirements: 



1- Njct > 2, with p T (j 2 ) > 40 GeV 

2. Ef iss > 30 GeV 

3. N\ cp = 2, where lep = e/fj, (isolated) and pr(h) > 20 GeV 

4. At least two jets tagged as b with pt > 50 GeV 

5. Only one of the two possible sets of pairings of the two leptons with the two leading 
6-jets should generate invariant mass values which are both less than 175 GeV. This 
cut was intended to reduce the experimental combinatorial background and was used 
only when constructing observables which required the pairing of leptons and jets 
from the same decay chain. 

Approximately 16100 (8300) events passed cuts 1-4 (1-5) respectively. Of these 15200 
(7400) were indeed events in which both hard W's decayed into a muon or electron. The 
remaining events contained at least one tau lepton decaying leptonically into e or /x. 

With the two bl pairs, each corresponding to the decay of a different top quark, one 
can construct the observables discussed in Section 5.1. These observables are m(b('' ,£(''), 



M C T{b,b'), McT(iJ') and M CT ([&4 [b'l']). Neglecting the mass of the 6-quark, the end- 
points in the distributions of these quantities lie at (from Eqns. (|5.2| )-( |5~5|) ): 



m-ft l) = vggEgfflfflZgwj _ 1SM GeV , m 

m(W) 

MEF(b,b>) = ^ ~^ {W) ee 135.0 GeV, (5.9) 
m[t) 

M max M = m\W)-m\u) ^ 8Q 4 GeV; 
m(W) 

MSrm , m = WHt) -^ {W) + m( t ) ( m ' (W J^ Hv) ) - 307.5 GeV.(5.11) 

Accounting for m(b) ^ translates into shifts of less than 0.1% in the end-point positions. 

We show in Figure |8| the distributions of the observables at parton- level and at detector- 
level for events passing the selection cuts in which both W 7 s decay into electron and muons. 
All contransverse mass observables have been corrected for transverse boosts according 



to the procedure discussed in Section 2.2. It can be seen that the end-point structures 
at parton-level are conserved at detector-level, modulo some smearing. The enhancement 
observed at the lower limit of the MctQ*, b') distribution is generated by the Jacobian peak 
at M C t = = 2m(b) discussed in Section | together with the M C t = M Cy effect of 

the boost correction procedure discussed in the same ssection. The Mqt distribution 
in Figure |8] is relatively unaffected by the Jacobian enhancement because M™™ = for 
massless leptons. The dilepton systems in these events receive large boosts from the bb' 
recoil however and so the boost correction procedure generates a prominent Mqt = Mq v 
peak at Mqt = 0. 
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Figure 8: Distributions of m(&W,^')) (top-left), Mct([M}, [b'i'}) (top-right), M C T(b,b') (bottom- 
left) and Mct{£,£') (bottom-right) for tt events passing the selection cuts where both leptons are 
generated directly from a W decay. The histograms show the parton-level distributions while the 
points with error-bars show the distributions after detector-level smearing. The expected end-point 
positions are indicated with vertical lines. The small populations of parton-level events lying beyond 
the expected end-points arise from the natural width of the W. 



In order to explore in an approximate manner the potential precision of mass measure- 
ments obtained with this technique, we fit the end-points of the distributions with a linear 
function smeared by detector resolution effects. Following Ref. 2S] we use a function f(x) 
given by: 

f(x) = -±=- J exp ( - )max{A(x / - x EP ), 0} dx 1 + a + bx. (5.12) 

Here x is the observable under consideration, x EP represents the end-point position, a 
represents the resolution of the assumed gaussian detector smearing, A is the slope of the 
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m(b,l) (GeV) M CT (b,b') (GeV) 

Figure 9: Distributions of m(b^'\ (left) and McT(b,b') (right) at detector-level for ti events 
passing the selection cuts. The dark (red) histogram indicates the distribution of events passing 
the selection where one of the two leptons is not directly produced in the decay of a W. The fit to 
the end-point function described in the text is shown. 




25 50 75 100 25 50 75 100 



M CT (1,1') (GeV) M CT (1,1') (GeV) 

Figure 10: Distributions of Mct(£,£') for A^-^m > (left), and A x n a b) < (right) at detector- 
level for ti events passing the selection cuts. The dark (red) histogram indicates the distribution of 
events passing the selection where one of the two leptons is not directly produced in the decay of a 
W. The fit to the end-point function described in the text is shown. 

distribution before smearing, and a and b are parameters describing an assumed linear 
background distribution. The latter distribution helps to take into account the effects 
of both combinatorial background from incorrect assignment of visible particles to decay 
chains and non-gaussian tails in the experimental resolution. 

For the observables m{b^\(S^) and Mct(&, 6'), the fits to the distributions for all 
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detector-level events passing the cuts are shown in Figure [|, with the irreducible background 
of ti events where at least one of the leptons is not directly produced in a W decay shown 
in grey (red). The fit function reproduces well the observed shape, and the value of the 
resolution parameter a obtained from the fit is in good agreement with the actual value of 
the smearing used in the detector simulation. 

The situation is somewhat more complex for the Mct(£,£') observable. In this case 
one has two populations. If -A^iab) > only a very small transverse boost correction 
is applied to Mct, using E^f = E cm . Therefore the experimental end-point resolution 
is to a good approximation just the resolution of the lepton pt measurement (of order 
1 GeV), plus the end-point smearing arising from the W natural width (of order 2-3 GeV). 
If A,. 

(lab) < however, Mct is corrected using the px of the hadronic recoil, resulting in a 
significantly larger resolution of order 9-10 GeV. The two configurations must therefore be 
fitted separately. The distributions are shown in Figure 1C, for A^ab) > (Ar(iab) < 0) on 
the left (right). From the figure one can also observe that the gradient of the A^iat,) > 
distribution near the end-point is smaller, and it was necessary in this case to fix the 
experimental resolution to 3 GeV in order to obtain an acceptable fit. 

Measurement of the end-point in the Mct ([&£], \b' £']) distribution presents further 




Figure 11: Two-dimensional distribution in the Mcr{\bP\, [b't 1 ]) versus m^ ax plane of detector- 
level ti events passing the selection cuts. The extremal values of the two observables, given by 
Eqns. ( p".llj ) and ( |5.8| ) are denoted by horizontal and vertical lines respectively. The dependence of 
Mgji* ([&£], [b'£' ]) on m^ ax given by Eqn. < $L7\) is denoted by the diagonal line. 
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Figure 12: Distributions of McT([b(], [b'l'}) for detector- level tt events passing the selection cuts 
after requiring additionally that m 2 nax < 20000 GeV 2 (left) and < 17500 GeV 2 (right). The 

vertical lines indicate the expected end-point positions. 



challenges due to the concave shape of the distribution near to the end-point, seen in 
Figure ^(top-right). For this end-point the assumption of a linear shape breaks down, pri- 
marily due to the depopulation of the Mcr([bl}, [b'l'}) versus m max plane near 



"max 

\2 



(m max (6, I)) seen in Figure 11. An alternative strategy for constraining the masses 



with M C T([b£], [b'l'}) would be to measure the dependence of Mffi([bt], [b'l'}) on m max , 
as discussed in Section 4T. This could be accomplished in practice by constructing 
Mct([W], [b'l'}) histograms of those events which pass a cut on m max . Unfortunately 
however this procedure is also complicated by the presence of concave end-points, as can 
be seen in Figure 

U. Further progress with this specific element of the contransverse mass 
technique will likely require a dedicated study of end-point shapes, which is outside of the 
scope of this paper. Because of these considerations we will not use the measurements of 
MjSji* ([&£], [b'l'}) in the following mass measurement study. Nevertheless such constraints 
could be useful for validating mass measurements obtained from the other observables. 
Similar considerations apply when attempting to use the non-boost-corrected Mcr(6, &') 
and Mct(£,£') observables to measure the dependence of M^ x (b,b') and M^ x (l,l') on 
Pb, as shown in Figure 13, and this technique will also not be considered further here. 

The results of the end-point fits are listed in Table |l|, where the first uncertainty is 
the statistical uncertainty from the MINUIT [30] fitting program for the chosen fit interval, 
the second is the systematic uncertainty obtained by varying the fit interval and the third 
uncertainty is the correlated systematic uncertainty derived from assumed energy scale 
uncertainties of 1% for 6-jets and 0.1% for leptons f3l[| . The quoted uncertainties should 
be considered approximate and could be improved with the use of better end-point fitting 
functions, for instance templates derived from Monte Carlo simulation studies. 

Based on the end-point measurement uncertainties listed in Table || it is possible to 
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Figure 13: Two-dimensional distributions in the Mcr(b, b') versus pb plane (left) and Mct{£,£') 
versus pb plane (right) of detector-level tt events passing the selection cuts. The dependence of 
M™^ x on pb given by Eqn. (4.9) is denoted by the curved line in each case. The small population of 
events lying beyond M™^ in the left-hand figure arises from the finite detector-level b-jet energy 
resolution, which degrades at lower energy. 



End-point 


Truth (GeV) 


Measured (GeV) 


m max (6,£) 






152.6 


152.8 ±1.7±1± 0.8 


Mffi x (b,b r ) 






135.0 


137.7±3.6±3± 1.4 


Mffi x (£,t') 


(Ar(lab) 


<0) 


80.4 


80.2 ± 0.5 ± 1 ± 0.1 


M^(£,£') 




>0) 


80.4 


81.2 ± 1.7±2±0.4 



Table 1: End-point positions in GeV. The first uncertainty is statistical, while the second and 
third are respectively the uncorrelated systematic and correlated energy scale uncertainties. The 
expected end-point positions from Eqns. (5.8), (5.9) and (5.10) are listed in the column labelled 
'Truth'. The assumed integrated luminosity is 3 fb _1 . 



evaluate the achievable precisions for measuring the masses of the top quark, W and neu- 
trino. We use the technique described in Ref. [p^| , where for each end-point measurement 
we generate a set of pseudo-experiments by sampling from a gaussian distribution centred 
on the nominal value of the end-point position of width equal to the estimated measure- 
ment precision. We assume that the measurements are uncorrelated with the exception of 
the energy scale uncertainties, which are assumed to be fully correlated. For each pseudo- 
experiment we calculate the value of m(t) according to Eqn. fl5.6| ) and hence calculate 
m(W) and miy) from Eqns. (5.9) and Q5.10 ). 

The distributions of the measured m(t) values and m(t) — m(W) mass differences are 
shown in Figure 14 for a set of 100000 pseudo-experiments. The precision of the top quark 
mass measurement is ~8 GeV, while the uncertainty on the measurement of m(t) — m(W) 
is ~2 GeV. A 95%(68%) upper limit on the neutrino mass of 30(16) GeV is obtained. These 
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Figure 14: Distributions of the calculated top mass (left) and m(t) — m{W) mass difference for 
100k experiments. The assumed statistics is 3 fb _1 . 



results may appear to be disappointing when compared with the ~1 GeV m(t) precision 
expected to be obtained at the LHC from semileptonic ti events for the same assumptions 
on 6-jet energy scale uncertainty Here however we have made no assumptions about 
the masses of the W or neutrino and so 8 GeV is the stand-alone precision with can 
be obtained with the technique. The end-point measurements used in this technique are 
primarily sensitive to mass differences and so if the mass of the W were assumed to be 
known the precision of the measurement of m(t) would improve to ~2 GeV, dominated 
by the systematics associated with the very crude end-point fitting function used for this 
study. 



5.3 A SUSY example 

Having demonstrated the proposed mass measurement technique with ti events let us now 
apply the same technique to a SUSY model generating events with a similar final state. 
An example of such a SUSY model is an MSSM model with a left-handed slepton doublet 
lighter than the chargino. In this case events with the decay chain 

q^QXi^qtv^qivXi (5-13) 

appearing in both legs of the event can be produced. The invisible sparticle at the end 
of the chain is in this case the sneutrino, since both of its decay products are undetected. 
This decay chain maps onto Eqn. ( p\l| ) with 5 = q, j3 = Xii<x=v, P = q and Q = I. 

The decay ( [5.13 ) is however not the only decay yielding the final state of interest, with 



a quark, a lepton and one or more invisible particles on each leg. We consider the following 



- 21 - 



decay chains: 



q ->■ qxf ->■ 9^ -> 9^X? (5-14) 
9 -> 9xf 9X?^ -> 9Xl^ (5-15) 
9 -> qxf -> 9X?^ -> 9Xi^ (5-16) 



For all of these decays, formulas (5.2)— (|5Tq) are valid, provided that m{a) in ( |5.2| ), 
and ( |5.5| ) is replaced by m m i n (a), defined as the minimum mass of the "pseudo-particle" 
composed of all of the invisible particles in the decay. Analytical expressions for ?Timin(^) 
in terms of the masses of the particles involved in the decays are given in Appendix B. 



These can be used in Eqns. (|5.2[) and (5.5). 



The case of Eqn. ( |5.4D deserves special comment. In this case for decay chains ( 5.14| ) 



(pTISD and ( gig ) one invisible particle is upstream of the lepton and one downstream. It 



is therefore not possible to correct for the upstream momentum, since it is not possible to 
separate the momentum of the neutrino and of the Xi- However, if one performs the boost 
correction assuming that all of the observed missing transverse momentum is downstream 
of the lepton, the distributions for M^ X (Q,Q') still possess end-points at positions given 
by Eqn. ([TJ) with the 

Wmmifi 1 ) values defined in Appendix B. 

Based on measurements of the end-point positions k%, and £3, one can calculate the 
masses of the squark and of the chargino independently from the decay mode of the X\ ■ 
This is a remarkable achievement, as it shows that it is possible to perform an absolute 
measurement of the chargino mass through its leptonic decay notwithstanding the fact that 
two invisible sparticles are present in the decay of each chargino. The interpretation of the 
meaning of the measured m(a) is dependent on the chargino decay mode, and the analysis 
based on the measurement of fei, &2> ^3 an d does not allow the discrimination of the 
different expressions given in Appendix B. 

It should be noted that the mass hierarchy and coupling structure implied by the 
presence of the decay chain given by Eqn. ( 5,13| ) imply also the existence, with a significant 



branching ratio, of the chain: 

9 -> 9X2 -> qd L -> 9^Xi, (5-17) 

which can be easily selected by requiring the presence of two leptons with the same flavour 
and opposite sign. This 'golden channel' for SUSY mass measurement at the LHC can 
potentially provide additional information regarding the masses of sparticles involved in 
the chargino decay chain. In the following we shall not show a complete analysis along these 
lines, which has been already developed in detail in e.g. Refs. |2fJ, 21, [3|]. We shall limit 



ourselves instead to showing for our example model that the invariant mass of two OS-SF 
leptons does give a characteristic end-point structure from the x\ ~~ & decay. Starting from 
this end-point, and combining it with the hard jets in the events it is possible to measure the 
masses of q, \\ & an d X? H0> H3> HI • ^ ^ s then straightforward to insert the measured values 
of these masses into Eqns. (^!^)-( |5.5|) and observe that under the assumption of decays 
( [BHD , flOp , and (|5l6|) the mass measurements from the two analyses are inconsistent. 
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Parameter 


Value (GeV) 


Parameter 


Value (GeV) 


771(5) 


520.0 


m(u L ) 


503.4 


m(e L ) 


157.1 


m(u e ) 


135.7 


m(xf) 


231.5 




232.0 


rra(xi) 


117.2 







Table 2: Masses of the relevant sparticles for the example MSSM point. 



End-point 


Position: Chain (ft. 13 


) 


Position: Chain ( 


5.14 


m max (q,£) 
M^(q,q') 

M^([ q £][q'£']) 


362.1 
396.9 
151.9 
727.2 


297.7 
396.9 
69.7 
668.5 



Table 3: Expected end-point positions in GeV for the decay chains (5.13) and (5.14) 



In order to explore the feasibility of this measurement technique, we used HERWIG 



6.5 [33, to generate events from a toy MSSM model incorporating the mass hierarchy 
present in decay chain ( |5.13j ). The masses of all of the squarks were set to 500 GeV and 
those of all the sleptons to 150 GeV. The three gaugino masses M±, M2 and M3 were 
set respectively to 120, 250 and 520 GeV while the higgsino mass parameter \i was set to 
400 GeV, tan (3 to 10, and rriA to 400 GeV. The trilinear couplings were set to zero. The 
relevant sparticle masses, as calculated by ISASUSY 7.75 |35| are listed in Table ^. The 
expected end-point positions for the chains ( |5,13j ) and ( |5.14| ) are listed in Table [3|, using 
the results of the discussion above regarding the treatment of m(a) in decay chain ( |5.14[ ). 

A total of 800 K events were generated, corresponding to an integrated luminosity of 
~12 fb _1 . The generated events were passed through the same parameterised detector 



simulation as for the top sample described in Section |5.2| . Events were selected with the 
following requirements: 

1. iV jet > 2, with p T (j x ) > 100 GeV and p T (h) > 50 GeV. 

2. Eip iss > 100 GeV. 

3. N\ ep = 2, where lep = e/fj, (isolated) and prih) > 20 GeV. The two leptons were 
required to possess different flavours. 

4. Veto all events with jets with pr(ji) > 20 GeV labelled as a 6-jet or r-jet. 

The veto on b and r labelled jets was applied to reduce the SUSY background from events 
containing top quark or r lepton decays. 

The requirement of leptons with different flavour reduces the signal by a factor two, but 
it is necessary, as the same-flavour signal is dominated by SUSY background events in which 



the two leptons are produced in the decay chain ( 5.17 ). This is demonstrated in Figure 15 
which shows the lepton-lepton invariant mass distribution for opposite-sign same-flavour 
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lepton pairs. The full line is the inclusive distribution, and exhibits the characteristic end- 
point structure from decay Q5.17 ) [20|; the grey (red) area indicates lepton pairs from the 
decays ( gig ) and (|T4]). Fi gure O demonstrates also that the characteristic lepton-lepton 
invariant mass end-point from the X2 ~~ ^ decay chain is indeed observable for this model. 

Following application of the selection cuts described above the only significant remain- 
ing background was from ti production and only this background is considered in the 
following. Approximately 9700 SUSY events passed cuts 1-4. Of these 7200 were indeed 
events in which both the muon and the electron were produced directly in the decay of 
a sparticle from the chains ( 5.13| ) or ( |5.14 ). In the remaining events at least one of the 
leptons was generated by the decay of a tau lepton produced in one of the two legs of the 
event. The number of ti background events was approximately 1400. 

The parton-level distributions for the observables m(q('\£^), Mct([q(], [q'£'])i Mct(q, q') 



and Mct{V-A') are shown in Figure 1£ for all events passing the selection cuts in which 
both legs in the event contain the chain (|5.13| ) or the chain ( pM4|) , The contransverse 
mass observables have been corrected for transverse boosts according to the procedure 
described in Section 12.21. For the reasons discussed in Section B.2I we shall not measure 



or exploit the Mcr([g€], [q'i']) end-points in the following analysis, nor the non-boost- 
corrected McT(b,b') and Mct(£,£') versus pb end-points. Nevertheless such constraints 
could be useful for validating mass measurements obtained from the other observables. 

The first step in the detector-level analysis is the calculation of the invariant mass 




100 

m(ll) (GeV) 



200 



Figure 15: Detector-level lepton-lepton invariant mass for opposite-sign same-flavour lepton pairs 
after the selection cuts except with the requirement of different flavours for the leptons. The full 
line is the inclusive distribution while the grey (red) area indicates the distribution for lepton pairs 
from the decays ( 5.13| ) and (5.17). 
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O iso 




O 100 — 




200 400 600 

M CT ([ql],[qT]) (GeV) 




> 
U 



200 400 
M CT (q,q') (GeV) 




100 200 

M rT (l,l') (GeV) 



Figure 16: Parton- level distributions of m{q^'\l^) (top-left), M C T{[qi), [q'i'}) (top-right), 
McT(q,q') (bottom-left) and Mct{(-,C) (bottom-right) for SUSY events passing the selection cuts 
where both leptons are generated by decay chain ( [5. 13 ). 



of each lepton with each of the two leading jets in the event. The distribution of the 



minimum of these two masses for each lepton is plotted in Figure 17 and displays an end- 
point at around 360 GeV, as expected from chain ( |5.13 ), The detector-level distributions 
of Mct(q,q') and Mct(£,£') are plotted in Figure |l^ and also display end-points at the 
positions expected for chain (|5 . 1 3 ) . Only the distribution of Mct(£,£') values for events 



with Ar(i a b) < is shown. This is due to the fact that the distribution at truth level for 
Ar(iab) < hits the nominal end-point, whereas that for A^ab) > runs out of statistics 
approximately 10 GeV below the nominal position (see Figure [Deleft)), leading to a biased 
fitted end-point position at detector level (see Figure |l9|(right)). This arises because only 
a very small boost correction (with Efg = E cm ) can be applied in the A x ^ ah ^ > case 
and hence the resulting corrected Mqt value is more conservative than in the A^iab) < 
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200 

m(q,l) (GeV) 



300 350 

m(q,l) (GeV) 



Figure 17: Detector- level distribution of the minimum value of m(gW,f (/) ) for SUSY events 
passing the selection cuts. The dark grey (blue) histogram indicates the tt background. The 
complete distribution is shown on the left while the region near the end-point is expanded on the 



right. The fit to the end-point function given in Eqn. (5.12) is shown 



End-point 


Truth (GeV) 


Measured (GeV) 


m max {q,£) 


362.1 


369.2 ±2±5± 1.5 




396.9 


401.7 ±4.8 ±5 ±4 


M%F{l,e) (A x <0) 


151.9 


149.3 ± 1.5 ±3± 0.8 



Table 4: End-point positions in GeV. The first uncertainty is statistical, while the second and 
third are respectively the uncorrelated systematic and correlated energy scale uncertainties. The 
expected end-point positions from Eqns. ( |5.2[ ), ( [j~3| ) and (5.4) are listed in the column labelled 
'Truth'. The assumed integrated luminosity is 12 fb _1 



case. This effect is most evident when consider ingMcr (A O because of the potentially 
large boosts generated by the recoiling W system. The same effect, although numerically 
less evident, is present also in the top quark analysis but is masked by the smearing of the 
end-point due to the W natural width. We choose here to use the larger of the two fitted 
end-point positions, which must be nearer to the true value. 

In order to explore the approximate potential precision of mass measurements obtained 
with this technique, we fit the end-points of the distributions with the smeared linear 
function given in Eqn. ( |5.12| ). The caveats associated with this technique discussed in 



Section 5.2 are also relevant here. In the case of the Mqt distribution we only fit 



the distribution of events with Aj,(i ab ) < 0, as discussed above. The fitted distributions are 



shown in Figures |17| and |18|. The results of the end-point fits are listed in Table [|, where 



the first uncertainty is the statistical uncertainty from the MINUIT [ 30 1 fitting program 
for the chosen fit interval, the second is the systematic uncertainty obtained by varying 
the fit interval and the third uncertainty is the correlated systematic uncertainty derived 
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200 400 
M CT (q,q') (GeV) 




50 100 

M CT (1,1') (GeV) 



Figure 18: Detector-level distributions of Mcr(q, q') (left) and Mct(£, £') with A x ( lab ) < (right) 
for SUSY events passing the selection cuts. The light grey (red) histogram is the distribution of 
events in which at least one of the two leptons was not produced directly from the decay of a 
sparticle. The dark grey (blue) area indicates the tt background. The fits to the end-point function 



given in Eqn. (5.12) are shown 



100 

M CT (1,1') (GeV) 



150 




50 100 150 

M CT (1,1') (GeV) 



Figure 19: Distributions of Mct(£, f) at parton-level (left) and at detector-level (right) for SUSY 
events passing the selection cuts. In the left-hand figure the dashed histogram is the distribution of 
events with A x (\ a \^ < while the full histogram is that of events with A^ab) > 0. In the right-hand 
figure all events possess A x ^ a ^ > and the shaded histograms and fitted curve are as for Figure |lj. 
The fitted end-point position is 133.6 GeV. 



from assumed energy scale uncertainties of 1% for jets and 0.1% for leptons As for 
the top study the quoted uncertainties should be considered approximate and could be 
improved with the use of better end-point fitting functions, for instance templates derived 
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from Monte Carlo simulation studies. 

If we use the measured end-point positions to calculate the masses of the sparticles 



using Eqns. ( |5.6| ), (5.3) and ( |5.4|) we obtain an uncertainty of 20 GeV on the absolute squark 
mass, an uncertainty of 6 GeV on the difference between the squark mass and the masses of 
the other sparticles, and an uncertainty of 3 GeV on the chargino-sneutrino mass difference. 
We have thus shown with a toy SUSY model that it is possible to achieve a stand-alone 
measurement of sparticle masses using the contransverse mass technique applied to events 
containing two symmetric sequential two-step two-body decay chains. 



6. Conclusions 

In this paper we have extended the contransverse mass technique for measuring the masses 
of pair-produced semi-invisibly decaying heavy particles so that it can be applied to events 
with non- negligible boosts of the CoM frame of the heavy states in the laboratory transverse 
plane. We have demonstrated the modified technique with case studies measuring the 
masses of the top quark, W and neutrino in fully leptonic ti events, and the masses of 
sparticles in SUSY events with a similar final state. The case studies presented here are in 
many respects more detailed than previous studies of alternative strategies and illustrate 
well the potential utility of the contransverse mass technique. 
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A. The connection between M T2 (x) and M C t 



Eqn. ( flip can be used to study the links between Mt2{x) an d Mct, as we shall now 
illustrate. In the process we shall obtain an approximate analytical expression for a boost- 
corrected version of M^ix)- The link between the two quantities in the absence of co-linear 
transverse boosts and for massless visible states was first discussed in Ref. [|36[j . 

First, observe that if we assume a value for m(a) and know Mqj? then we can solve 
Eqn. (4.1) for m{5). For any given event we do not know M™^ x however, but rather Mqt- 
Let us therefore substitute Mqt for Mq^ to obtain the following solution 



m so i n (5) = (x 2 + A T 



1 + 



4x 2 



2At — m 2 {v\) — m 2 {v2) 



m 2 {vi)m 2 {v2) 



(Al) 
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where At = \M CT — m 2 (vi) — m 2 (v2)]/2 and x is the assumed value of m(a). This is 
identical to the analytical expression for the 'balanced' solution of Mt2{x) i n the absence 
of co-linear transverse boosts, which was identified in Refs. |37| j7|. In general m so i n (5) 
need not be bounded by m(5), because At and hence Mct appears in the denominator 
inside the square-root causing m so \ n (5) not to be a monotonically increasing function of 
Mct- Therefore a value of Mct < Mqj? need not generate a value of m so \ n (8) which 
is less than the value obtained with Mct = M™^ x . Note that if x = then m so \ n (5) 
is nevertheless bounded by m{$) because in this special (5) is a monotonically 

increasing function of Mct- 

In general m so i n (S) is not invariant under co-linear transverse boosts, because it de- 
pends on Mct which is also not invariant. One might consider therefore whether it is 
possible to correct m so \ n (5) for such boosts in a similar manner to that used to correct 



Mct i n Section 2.2. The Mct boost-correction procedure minimises Mct with respect to 
the possible boosts, however this does not necessarily minimise m so \ n (S) because it is not 
in general a monotonically increasing function of Mct- Boost-correction can be performed 
however, as we know not only the minimum possible value of Mct in the <5i<52 CoM frame 
but also the maximum possible value, which is given by the maximum of the Mct values 
obtained by boosting v\ and V2 with respectively (3 = Pb/ E cm and f3 = Pb/E (see Fig. |||) . 
We therefore know the range of CoM frame Mct values which could have occurred in the 
event. 



To proceed we find the turning-points of Eqn. (Al), which lie at the following four 



values of M CT : 

2 _ x[3m 2 (tQ + m 2 (^)] ± 2m(vi)[m 2 (vi) + m 2 (v j)] 

IV1 CT — — — , {AZ) 

X ± m(Vi) 

where i,j = 1,2 (i ^ j). We then take the minimum of 3-11 the ^soln 

(5) values obtained at 

turning-points where Mct hes within the allowed range identified above. We finally take 
the minimum of this m so i n (<5) value and those obtained with Mct set to its extrema. The 
result is a boost-corrected value of m so i n (8). 

Now although m so i n (5) is not in general bounded by m(5) the full analytical expres- 
sion for Mt2{x) h 1 the absence of co-linear transverse boosts [^, 0] is so bounded. This 
is because in those cases where m so i n (5) > m(S) Mt2(x) takes on an alternative value 
corresponding to an 'unbalanced' solution (see e.g. Eqn. (54) in Ref. [p7| ). The quantities 
upon which the decision to switch to such an alternative value rests are not themselves 
boost-invariant. Development of an appropriate boost-correction procedure for these quan- 
tities consistent with the parallel boost-correction of m so \ n (5) requires more work. In the 
absence of such a correction procedure we can choose to neglect the boost when making 
this decision, to obtain an approximate analytical form for a boost-corrected version of 

M T 2(x)- 

Note that the quantity calculated here is not the same quantity as the conventional 
Mt2(x) used elsewhere, which is boost-independent but currently does not possess a general 
analytical form. The boost-correction procedure, even if exact (i.e in the absence of the 
approximation mentioned in the previous paragraph), leads to a quantity which is not in 
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general equal to the conventional Mt2(x)i although equality is obtained in the absence of 
co- linear transverse boosts. 

The procedure described above for calculating this 'boost-corrected M^ixY quan- 
tity is implemented in the boost-correction package described in Section |2.2| , available at 
http : //projects .hepforge . org/mctlib. 

B. Invisible pseudo-particle masses for chargino decay chains 

If a heavy sparticle decays via a chain which produces multiple invisible final state particles 
then the values of m max (P,Q) and M™| X ([PQ], [P'Q 1 ]) can be calculated by constructing 
an aggregate 'pseudo-particle' a from the invisible particles. The minimum value of m(a), 
^mm(aO) can then be used in end-point formulae, as described in Section 

If a chargino decays through xt ~^ Xi^ ± ~~ > Xi^i^ then a = [Xi^] and m mm (a) is 
given by: 

m 2 min (a) = m 2 (x°i) + m(W) [E(&) - p] J (Bl) 



E(W)+p' 



where 



\j [m 2 (xf ) - rn 2 {x\) ~ m 2 {W)] 2 - Am 2 (xi)m 2 {W) 
2m(xf) 



P = 1 ^-7^ > (B2) 



jp fV( r\ rn 2 {xf)-m 2 (xi) + m 2 (W) 

E{W) = ^ixf) ' (B3) 



w~ck = m 2 (xf) + m 2 (x°i) ~ m 2 (W) 
2m(xf^ 



and we have assumed m{y) = 0. If however the chargino decays through xt ~^ v^ 1 ~^ 
vl~^x\i as in decay chain fl5,14|) , then: 

m min (a) = -~- , (B5) 

while if the chargino decays through xt — > t^v — > i^uXx, as in decay chain ( |5.13|) , then 
a = v and m m i n (a) is fixed to m(a) given by: 

m(a) = m(u). (B6) 

Finally, if the chargino decays through the three-body decay xt ~~ ^ X?^ ±I/ then a = [Xx l/ ] 
and m m i n (a) is given by: 

m min (a) = m(xi). (B7) 
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